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Abstract. We discuss the possibility that microquasar jets may be powerful emitters 
of TeV neutrinos. We estimate the neutrino fluxes produced by photopion production 
in the jets of a sample of identified microquasars and microquasar candidates, for which 
available data enables rough determination of the jet parameters. We demonstrate that 
in several of the sources considered, the neutrino flux at Earth, produced in events similar 
to those observed, can exceed the detection threshold of a km 2 neutrino detector. Sources 
with bulk Lorentz factors larger than those characteristic of the sample considered here, 
directed along our line of sight may be very difficult to resolve at radio wavelengths and 
hence may be difficult to identify as microqusar candidates. However these sources can 
be identified through their neutrino and gamma-ray emission. 



1. Introduction 

The composition of microquasars jets is yet an open issue. The synchrotron emis- 
sion both in the radio and in the IR is consistent with near equipartition between 
electrons and magnetic field, which is also implied by minimum energy consider- 
ations [p. However, the dominant energy carrier in the jet is presently unknown 
(with the exception of the jet in SS433). A possible diagnostic of hadronic jets 
is emission of TeV neutrinos ||. As shown in reference 0, for typical micro- 
quasar jet parameters, protons may be accelerated in the jet to energies in excess 
of ~ 10 16 eV. The interaction of these protons with synchrotron photons emit- 
ted by thermal electrons is expected to lead to 1-100 TeV neutrino emission. The 
predicted fluxes are detectable by large, km 2 -scale effective area, high-energy neu- 
trino telescopes, such as the operating south pole detector AMANDA and its 
planned 1 km 2 extension IceCube ||, or the Mediterranean sea detectors under 
construction (ANTARES §; NESTOR @) and planning (NEMO @ and § for 
a recent review). 

In this paper we consider a class of identified Galactic microquasars with ei- 
ther persistent jets or documented outbursts. For each source we provide, for 
illustrative purposes, our model prediction for the neutrino flux that should have 
been emitted during particular events, using radio data available in the literature. 
Although the temporal behavior of many of these sources may be unpredictable, 
we demonstrate that some of the sources could have been detected by a neutrino 
telescope with effective area larger than km 2 (in some cases even 0.1 km 2 ) had 
such a detector been in operation during the time of the recorded events and, 
therefore, propose that they should be potential targets for the planned neutrino 
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telescopes. In addition we consider a list of XRBs thus far unresolved at radio 
wavelengths, that are believed to be microquasar candidates. In §|2] we briefly 
discuss the neutrino production mechanism in microquasars. In §|| we use obser- 
vational data available for each source to estimate the jet parameters, and then 
use these parameters to derive the expected neutrino flux. The number of neu- 
trino induced muon events in km 2 -scale neutrino telescopes is derived in §|I[ The 
implications of our results are briefly discussed in §g. 

2. Internal shock model for microquasars 

In this section we give a brief outline of the model proposed by Levinson &; 
Waxman (2001) Q for production of neutrinos in microquasars with the purpose 
to introduce the parameters relevant for the present analysis. The model assumes 
that on sufficiently small scales (< 10 11 cm), unsteadiness of the jet leads to 
the formation of internal shocks that can accelerate protons and electrons to a 
power law distribution. For typical jet parameters, the maximum proton energy is 
roughly 10 16 eV in the jet frame. Protons can interact with both the external X- 
ray photons emitted by the accretion disc, and the synchrotron photons produced 
inside the jet by the shock heated electrons, leading to pion production and 
consequent neutrino emission. In order for photomeson production to take place, 
the comoving proton energy must exceed the threshold energy for A-resonance, 
~ 10 eV for interaction with the external photons and ~ 10 13 eV for interaction 
with synchrotron photons. 

Charged pions produced in photo-meson interactions decay to produce neu- 
trinos, 7r^ — > fi + + — > e + + v e + + v^. In a single collision, a pion is 
created with an average energy that is ~ 20% of the proton energy. This energy 
is roughly evenly distributed between the final tt + decay products, yielding a v 
energy that is ~ 5% of the proton energy. The neutrino signal is dominated by 
neutrinos in the energy range of 1 to 100 TeV. The fraction of proton energy 
converted into pions, f n , depends on the jet Lorentz factor, T, and on the kinetic 
luminosity of the jet, Lj et (see Q for details). 

The expected neutrino flux (energy per unit area per unit time) of ^ 1 TeV 
muon neutrinos at Earth from a jet ejection event is p] 

where 5 = [T(l — /3cos#)] _1 is the Doppler factor of the jet (9 is the angle 
between the jet axis and our line of sight), rj p ~ 0.1 is the fraction of Lj e t carried 
by accelerated protons, and D is the source distance. 

3. Jet parameters and neutrino fluxes 

In order to determine the neutrino flux from (|lj) the main quantity to be estimated 
is the kinetic power of the jet. We have used two different methods for the sources 
with the relativistic jets resolved in the radio band (resolved microquasars) and 
for the XRBs observed as point like sources, such as GX 339-4, henceforth referred 



to as unresolved microquasars 10]. 
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3.1. Resolved microquasar s 

In events that have been monitored with sufficiently good resolution, it is often 
possible to obtain a rough estimate of the characteristic source parameters, in 
particular the bulk speed of the jet, the angle between the jet axis and the sight 
line, and the size of the emitting blob. Lj e t is estimated in the following way: 

for a flux density S u oc v^ 1 / 2 of the radio source, implying an electron energy 
distribution dn e /de e oc e -2 , we estimate the minimum energy carried by electrons 
and magnetic field, obtained for a magnetic field [e.g. 



£?* = 3.6 



hl(7max/7mi ^ ' 
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where v = ug ■ 10 9 Hz, I = lis ■ 10 15 cm is the size of the emission region, 7 n 
and 7 m i n are the maximum and minimum electron random Lorentz factors as 
measured in the jet frame, and Tg = 10 6 Tge K is the brightness temperature. 
The jet power then satisfy, 

L jet > 0.3c(r/£g 2 . (3) 

For the numerical estimates that follow we conservatively assume 7max/7min = 
100. In Table |l| we report our estimates of the jet power and of the neutrino 
flux expected at Earth, calculated using equations n] and || and adopting the 
source parameters quoted in literature. The two values of the periodic source 
LS I +61°303 (P~26.5 days) refer respectively to bursting and quiescent states 
observed by Massi et al. (2001) |TT| ]. The distance of V4641 Sgr, and therefore 
the jet parameters, is un certai n. W e report calculations for both the values of D 
present in literature ||12| | and |13|]. Several authors pointed out that the kinetic 
energy output of the Ss433 jets can influence the radio emission of W50, moreover 
SS433 is the only microquasar that shows a strong H a -line emission from the jets 
|j~4| , [15| , 16]. In our estimate we assume the conservative value of Lj et ~ 10 39 
erg/sec suggested by Margon (1984) [p7|. 



3.2. Unresolved microquasars 

For the sources whose jet has not been resolved, we cannot deduce the value of 
Lj e t from Eq- (|3j) , since T, 5 and the size of the jet are not known. We follow 
instead a different line of argument. We estimate the jet synchrotron luminosity, 
assuming a spectral index or ~ 0.5 for the emitted radiation [l(J: 

L syn = 4irD 2 1 S u u high , (4) 

1 — OtR 

where Vhigh is the highest observed frequency of synchrotron emission, and S Uhigh 
the flux density emitted at this frequency. Assuming that the total jet synchrotron 
luminosity is a fraction r\ r of the jet kinetic energy, carried by relativistic electrons 
and magnetic field, we estimate the total Ljet as: 

Ljet = Ve lr lr ^sj/m (5) 

where r\ e denotes the fraction of bulk energy converted internal energy of elec- 
trons. For a typical value of an ~ 0.5, that implies electrons do not cool fast on 
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scales that are resolved by the VLA. Assuming that emission at Vhigh originates 
from the same radius as the radio emission, for the typical parameters inferred 
for the resolved microquasars (B of tens of mG corresponding synchrotron fre- 
quency of about 10 GHz), electrons radiating at Vhigh ~ 10 14 cannot lose more 
than r) r ~ 0.1 of their energy. The energy neutrino flux expected at Earth is: 

_1 Lj e t/8 _ 1 rjp _x . . 

~ 2 UVp 4ttD2 ~ 16(1 - a R ) r,/^ bv ^hi g h- W 

We have neglected in eqs. (Q-f|) corrections associated with the relativistic ex- 
pansion of the jets. However, since these corrections are the same for both the 
synchrotron and neutrino emission, the estimate of neutrino flux in eqs. (|q) is 
independent of such corrections. In Table [j] we quote, for a sample of unresolved 
microquasar candidates, our estimates of Lj et and f v , calculated from Eq. (JsJ) 
and Eq.(|). 

Recent observations with the VLBA array have confirmed the nature of mi- 
croquasars for the source Cygnus X-l, revealing an extended jet-like feature ex- 
tending to ~ 15 mas from a core region, with an opening angle of < 2° |fL8| . The 
authors suggest for this source a bulk motion velocity of j3 ~ 0.75 and a viewing 
angle of 9 ~ 40°, corresponding to a Lorentz factor of V ~ 1.51 and a Doppler 
boost factor 5 ~ 1.55. Adopting these jet parameters and a quasi-steady emission 
of S u ~ 5 mJy at v = 8.4 GHz, the kinetic luminosity, calculated from Eq. (|3|), 
is Lj et ~ 3.53 • 10 36 erg/sec. The new estimate of the Lj et for Cygnus X-l is 
therefore consistent with the value in Tab. [I], obtained from the Eq. (||) for the 
unresolved sources. 



4. Muon events expected in a km 2 -scale detector 

The detection of TeV neutrino fluxes from microquasars could be the first achiev- 
able goal for proposed underwater (ice) neutrino telescopes. In this section we cal- 
culate the rate of neutrino-induced muon events expected in a detector with an 
effective area of 1 km 2 . Since the signal is expected to be dominated by neutrinos 
of energy E v > 1 TeV, for which the detection probability is P Ufl ~ l._3-10 -6 E v ^eV 



1 19], we estimate the number of neutrino-induced muon events as 



~ O^np^US^D^Lj^ssA^AU (7) 

where A km 2 is the effective detector area in km 2 units, and Atj the duration of 
the observed burst measured in days. It is seen from eq. [?| that in microblazars 
(microquasars having their jets pointing in our direction), even single outbursts 
may be easily detectable by km 2 detectors if the total energy release is of the order 
of that seen in GRS 1915 and GRO 1655 (L jet At ~ 10 43 ergs). In Table [l] we 
report the number of events expected in a detector with A e jj = 1 km 2 , during 
the bursts considered in section [31. In the case of persistent sources number of 
neutrino induced muon events is calculated for a 1 year period. In the same 
table, we also report the number of background atmospheric neutrino events 
collected in such a detector during the time At, assuming a neutrino spectrum 
4>v,bkg ~ 10 -7 E~^y i 'cm 2 sec sr for E v > 1 TeV and for a detector angular 
resolution of 0.3°. 
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5. Discussion 

There are large uncertainties involved in the derivation of the jet parameters for 
most of the sources listed in table 4. The best studied cases are perhaps GRS 
1915+105, GRO J1655-40, and SS433. Nonetheless, we have demonstrated that if 
the jets in microquasars are protonic, and if a fraction of a few percent of the jet 
energy is dissipated on sufficiently small scales, then emission of TeV neutrinos 
with fluxes in excess of detection limit of the forthcoming, km 2 scale, neutrino 
telescopes is anticipated. 

The present identification of microquasars, and the inferred distribution of 
their jet Lorentz factors, may be strongly influenced by selection effects. It is 
quite likely that the class of Galactic microquasars contains also sources with 
larger bulk Lorentz factors and smaller viewing angles, which should emit neu- 
trinos with fluxes considerably larger than the extended microquasars discussed 
in this paper. Such sources may be identified via their gamma-ray emission with, 
e.g. AGILE and GLAST, or by their neutrino emission. The gamma-rays should 
originate from larger scales where the pair production opacity is sufficiently re- 
duced. Predictions for AGILE and GLAST will be discussed in a forthcoming 
paper |2C]. There are currently about 280 known XRBs [21, 22], of which ~ 50 
are radio loud. These may also be potential targets for the planned neutrino de- 
tectors. Our results quoted in Table |l] are consistent with experimental upper 
limits on neutrino fluxes from point sources set by MACRO [23|| , AMANDA [24] 
and SuperKamiokande p5|. 
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TABLE 1. Estimates of the kinetic luminosity of the jet and neutrino flux for known microquasars. N M 
is the predicted number of muon events in a 1km 2 detector during the burts duration At and N^^kg 
is the expected number of background atmospheric neutrino events, assuming angular resolution of 
0.3°. 



Source name 
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CI Cam 
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Cygnus X-3 M 
LS 5039 [H 
GRO J1655-40 
GRS 1915+105 
Circinus X-l |2E 
LS I 61°303 
LS I 61°303 
XTE J1550-564 @ |7) 
V4641 Sgr [pi 
V4641 Sgr p , pi 
Scorpius X-l |3J, |k 
SS433 @ 
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41] 



43] 




3.62-10 37 
3.86-10 38 
1.45-10 36 
4.35-10 37 
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